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Abstract 

A search is presented for pair production of a new heavy quark (Q) that decays into a W 
boson and a light quark (q) in the final state where one W boson decays leptonically (to an 
electron or muon plus a neutrino) and the other W boson decays hadronically. The analysis 
is performed using an integrated luminosity of 20.3 fb“^ of pp collisions at = 8 TeV col¬ 
lected by the ATLAS detector at the LHC. No evidence of QQ production is observed. New 
chiral quarks with masses below 690 GeV are excluded at 95% confidence level, assuming 
BR(2 ^ Wq) = 1. Results are also interpreted in the context of vectorlike quark models, 
resulting in the limits on the mass of a vectorlike quark in the two-dimensional plane of 
BR(G ^ Wq) versus BR(2 ^ Hq). 
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1 Introduction 


The recent observation of the long-predicted Higgs boson by the ATLAS [1] and CMS [2] collaborations 
now completes the Standard Model (SM). In spite of its success, the SM cannot account for dark matter 
and the matter/antimatter asymmetry in the Universe and also fails to provide insight into the fermion 
mass spectrum, nonzero neutrino masses, why there are three generations of fermions, or why parity 
is not violated in the strong interaction. Furthermore, the observed Higgs boson is unnaturally light, 
requiring fine tuning to cancel radiative corrections that would naturally result in a mass many orders of 
magnitude larger, a discrepancy known as the hierarchy problem [3]. 

A variety of models have been proposed to address the various shortcomings of the SM. For example, 
a primary motivation for supersymmetry (SUSY) is to solve the hierarchy problem [4]. In SUSY mod¬ 
els, the quadratically divergent radiative corrections to the Higgs-boson mass due to SM particles are 
automatically canceled by the corrections from the supersymmetric partners. Models such as Little 
Higgs, Composite Higgs, and Topcolor take a different approach, proposing that electroweak symmetry 
breaking happens dynamically as the result of a new strong interaction [5-9]. Grand unified theories 
(GUTs) go further, unifying the electroweak and strong forces by proposing that the SM gauge symmetry 
SU (3)c X 5'f7(2)L X f7(l)Y is the low-energy limit of a single fundamental symmetry group such as SOCIO) 
or Fg [10, 11], which could potentially explain the observed spectrum of fermions and even provide in¬ 
sight into the unification of the electroweak and strong forces with gravity. A feature in many of these, 
and other models [12-14], is the prediction of vectorlike quarks (VLQs), hypothetical spin-^ particles 
that are triplets under the color gauge group and have identical transformation properties for both chiralit¬ 
ies under the electroweak symmetry group. Furthermore, massive VLQs would respect gauge invariance 
without coupling to the Higgs field. This allows VLQs to avoid constraints from Higgs-boson produc¬ 
tion [15]; if the Higgs sector is minimal, these constraints rule out additional chiral quarks. However, 
some two-Higgs-doublet models are able to avoid those constraints and accommodate a fourth generation 
of quarks [16]. 

In the models of interest, the VLQs have some mixing with the SM quarks, allowing them to decay to SM 
quarks and either a W, Z, or Higgs boson; however, the exact nature of the coupling depends on the model. 
For example, in composite Higgs models, the VLQs are involved in a seesaw mechanism with the SM 
quarks, so the lightest VLQ couples almost exclusively to the heaviest SM quarks (t- and (^-quarks) [6]. 
However, there are also models that predict TeV-scale VLQs that could preferentially decay to light SM 
quarks (q = u,d,s, or c) [10, 11, 17]. For example, the left-right mirror model (LRMM) [17] predicts 
three generations of heavy “mirror” quarks, with the lightest mirror generation coupling to the lightest SM 
generation. The two lightest mirror quarks could be pair-produced at the LHC via the strong interaction 
and would then decay to Wq, Zq, or Hq {q - u ox d). The LRMM would provide an explanation for 
tiny neutrino masses, parity violation in weak interactions, parity conservation in strong interactions, and 
could be the first step toward uncovering the symmetry structure of a GUT. Another model predicting 
VLQs that decay to light quarks is the E(, GUT with isosinglet quarks [10, 11]. In this model, after the 
Fg symmetry is broken down to the SM group structure, VLQ partners to the d-, s-, and fj-quarks are 
predicted. If the VLQs have the same mass ordering as their SM partners, the lightest VLQ would couple 
predominately to first-generation SM quarks {q - uox d). The values for the branching ratios to the three 
decay modes {Wq, Zq, Hq) depend on parameters in the model. The values in the Fg isosinglet model 
range from approximately (0.6, 0.3, 0.1) to (0.5, 0.25, 0.25), while the LRMM allows branching ratios 
from approximately (0.6, 0.4, 0) to (0, 0, 1), depending on the VLQ mass and mixing angles. 
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If such new quarks exist, they are expected to be produced predominantly in pairs via the strong in¬ 
teraction for masses up to 0(1 TeV) in LHC collisions with a center-of-mass energy of 8 TeV. Single 
production of a new heavy quark, Q, would be dominant for very high quark masses, but the production 
rate is model-dependent and could be suppressed if the coupling to SM quarks is small. To date, there 
have been two analyses of LHC data sensitive to VLQs that decay to light quarks, both using 1.04 fb“^ of 
ATLAS data ai ^[s -1 TeV: one search for pair-production [18] and one for single production [19]. The 
pair-production search set a lower limit on the Q mass of 350 GeV at 95% confidence level, assuming 
BR(2 —> Wq) = 1. Such a signal was also ruled out by the Tevatron for masses up to 340 GeV [20]. 



Figure 1: Leading-order Feynman diagrams for QQ —> WqWq production at the LHC. 


This paper presents a search for new heavy quarks that couple to light SM quarks using data collected 
by the ATLAS detector. The analysis focuses on model-independent pair production of a new heavy 
quark and its antiparticle, which then decay through a charged-current interaction to a final sfafe wifh a 
single elecfron or muon, missing fransverse momenfum and lighf SM quarks, making if complemenfary 
fo searches for VLQs fhaf decay fo fhird-generafion quarks [21-28]. The dominanf Feynman diagrams 
for fhe signal process are shown in Fig. 1. Background rejection is achieved fhrough even! topology 
and kinemafic requiremenfs. In particular, a kinemafic variable mofivafed by fhe spliffing scale of a 
heavy objecf info daughfer parficles [29, 30] is used fo reduce fhe background when selecting fhe pair of 
Jefs consisfenf wifh a hadronically decaying W boson. The sfalislical inferprefafion of fhe dafa uses fhe 
invarianf mass disfribufion of fhe Q candidafe formed from fhe hadronically decaying W boson and a lighf 
quark. The resulfs are inferprefed bofh in fhe confexf of VLQs and a chiral fourth-generafion quark. 


2 The ATLAS detector 

The ATLAS defecfor [31] af fhe LHC covers nearly fhe enfire solid angle around fhe inferacfion poinf 
(IP).^ If consisfs of an inner fracking defecfor surrounded by a fhin superconducfing solenoid, elecfro- 
magnefic and hadronic calorimefers, and a muon specfromefer incorporafing fhree large superconducting 
foroid magnefs. 

The inner fracking defecfor sysfem is immersed in a 2 T axial magnefic field and provides charged-parficle 
fracking in fhe range \q\ <2.5. A high-granularify silicon pixel defecfor covers fhe inferacfion region and 
fypically provides fhree measuremenfs per frack. If is followed by a silicon microsfrip fracker, which 
usually provides four fwo-dimensional measuremenf poinfs per frack. These silicon defecfors are com¬ 
plemented by a fransifion radiafion fracker, which enables radially exfended frack reconsfrucfion up fo 
\q\ - 2.0. The fransifion radiafion fracker also provides elecfron identification information based on 


* ATLAS uses a right-handed coordinate system with its origin at the nominal IP in the center of the detector and the z-axis 
along the beam pipe. The x-axis points from the IP to the center of the LHC ring, and the i/-axis points upwards. Cylindrical 
coordinates (r, (jy) are used in the transverse plane, 0 being the azimuthal angle around the z-axis. The pseudorapidity is defined 
in terms of the polar angle 6 as ?/ = - In tan(0/2). Angular distance is measured in units of hR = + (A^)L 
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the fraction of hits (typically 30 in total) above a higher threshold for energy deposits corresponding to 
transition radiation. 

The calorimeter system covers the pseudorapidity range \r]\ < 4.9. Within the region \q\ < 3.2, electromag¬ 
netic calorimetry is provided by barrel and endcap high-granularity lead/liquid-argon (LAr) calorimeters, 
with an additional thin LAr presampler covering It/I < 1.8, to correct for energy loss in material upstream 
of the calorimeters. Hadronic calorimetry is provided by the steel/scintillator-tile calorimeter, segmen¬ 
ted into three barrel structures within |7/| < 1.7, and two copper/LAr hadronic endcap calorimeters with 
\r]\ < 3.2. The solid angle coverage is completed with forward copper/LAr and tungsten/LAr calorimeter 
modules optimized for electromagnetic and hadronic measurements respectively. 

The muon spectrometer comprises separate trigger and high-precision tracking chambers measuring the 
deflection of muons in a magnetic field generated by superconducting air-core toroids. The precision 
chamber system covers the region \q\ < 2.7 with three layers of monitored drift tubes, complemented by 
cathode strip chambers in the forward region, where the background is highest. The muon trigger system 
covers the range \r]\ <2.4 with resistive plate chambers in the barrel, and thin-gap chambers in the endcap 
regions. 

A three-level trigger system [32] is used to select events for offline analysis. The Level-1 trigger is 
implemented in hardware and uses a subset of detector information to reduce the event rate to at most 75 
kHz. This is followed by two software-based trigger levels that together reduce the event rate to about 
400 Hz. 


3 Signal and background samples 

The pair-production cross section for a new heavy quark ranges from 12pb for a 300 GeV quark to 
21 fb for a 800 GeV quark. It was calculated at next-to-next-to-leading-order (NNLO) in QCD, including 
resummation of next-to-next-to-leading logarithmic (NNLL) soft gluon terms, with Top-t-t2.0 [33-38], 
using the MSTW2008 NNLO [39, 40] set of parton distribution functions (PDFs) and as = 0.117. The 
PDF and as uncertainties were calculated using the PDF4LHC prescription [41] with the MSTW2008 
68% CL NNLO, CTIO NNLO [42], and NNPDF2.3 5f FFN [43] PDF sets. The uncertainties in the 
prediction stem from scale variations and the PDF-i-as uncertainty and range from approximately 11% to 
12% for masses from 300 GeV to 800 GeV. 

VLQ signal samples were simulated with the tree-level event generator CompHEP v4.5. 1 [44] at the parton 
level with the CTEQ6L1 EO PDE set [45] and with the QCD scale set to the mass of the heavy quark, 
mg. The generated events were then passed into Pythia v8.165 [46, 47] for hadronization and parton 
showering. The VEQ signal samples were produced for values of mg ranging from 300 GeV to 800 GeV 
in 100 GeV steps. This range is motivated by the previous limit at 350 GeV and the expected sensitivity 
of this analysis to masses up to approximately 700 GeV. Although the analysis is targeting the Q —> Wq 
decay, there is also sensitivity to the neutral-current decays Q ^ Zq and Q Hq (e.g., QQ WqZq 
€vqqqq) and events were generated for all six decay combinations. In addition to the VEQ signal samples, 
a set of fourth-generation chiral-quark signal samples were generated with Pythia v8.1 using the MSTW 
2008 EO PDE set. The kinematics of the chiral-quark signal samples, which only contain QQ WqWq, 
are compatible with the VEQ samples when requiring BR(2 ^ Wq) = 1. Therefore, the more generic 
VEQ samples are used for the statistical analysis, with the sample corresponding to mg - 700 GeV and 
BR(2 ^ Wq) = 1 used to represent the signal in tables and figures, unless noted otherwise. 


4 


The background originates mainly from IT-boson production in association with jets, VT+jets, with lesser 
contributions from top quark pair production {tt), Z+jets, single-top, diboson, and multijet events. The 
VT-i-jets and Z-i-jets samples were produced using Sherpa vl.4.1 [48] with the CTIO PDF set, taking 
the c- and ^-quarks as massive. Samples of tt and single-top events were generated with Powheg-Box 
3.0 [49, 50] interfaced to Pythia v6.426 [47] using the Perugia2011C set of tunable parameters [51] 
for the underlying event and the CTEQ6L1 PDF set. Diboson production was modeled using Alpgen 
v2.13 [52] interfaced to Herwig v6.520 [53] with the CTEQ6L1 PDE set, except for the lepton-i-jets final 
state {WV —> ivqq with V = W,Z), which used the Sherpa vl.4.1 event generator with the CTIO PDE 
set. The contribution from multijet events originates from the misidentification of a jet or a photon as 
an electron, or from the semileptonic decay of a b- or c-quark, and the matrix method [54] is used to 
determine the kinematic distributions for the multijet background. 

The W/Z-i-jets and multijet backgrounds are normalized to data and a data-driven correction is applied 
to the transverse momentum p^ of the boson as described in Sec. 5.1. The tt cross section is determ¬ 
ined by the Top-i-i- prediction, with the top-quark mass taken to be 172.5 GeV. The single-top samples 
are normalized to the approximate NNEO theoretical cross sections [55-57] calculated using the MSTW 
2008 NNEO PDE set. The diboson background processes are normalized to NEO theoretical cross sec¬ 
tions [58]. 

All simulated samples include multiple pp interactions per bunch crossing and simulated events are 
weighted such that the distribution of the average number of interactions per bunch crossing agrees with 
data. The generated samples are processed through a simulation [59] of the detector geometry and re¬ 
sponse using Geant4 [60], then reconstructed using the same software as used for data. Simulated events 
are corrected so that the object identification efficiencies, energy scales, and energy resolutions match 
those determined in data control samples. 


4 Event selection 

The data analyzed in this search were collected with the ATEAS detector in 2012 and correspond to an 
integrated luminosity of 20.3 fb“^ Data quality requirements are applied to remove events with incom¬ 
plete, corrupted, or otherwise compromised subdetector information. Events are required to pass a single¬ 
electron or single-muon trigger. The pt thresholds are 24 GeV and 60 GeV for the electron triggers and 
24 GeV and 36 GeV for the muon triggers. The lower-threshold triggers include isolation requirements 
on the candidate leptons, resulting in inefficiencies af higher pj that are recovered by the higher-threshold 
triggers without an isolation requirement. 


4.1 Preselection 

The basic object selection is called preselection and requires exactly one charged lepton (electron or 
muon), at least four jets, and large missing transverse momentum as described below. The criteria 

are similar to those used in recent ATEAS top-quark studies [61], except that this analysis requires that 
there are no jets identified as originating from a fj-quark. The expected and observed numbers of events 
after preselection are shown in Table 1. There is negligible sensitivity to heavy quark production because 
the signal expectations for all masses are much smaller than the uncertainty on the background, dominated 
at this stage by systematic uncertainties. 
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Charged-lepton requirements 


Electron candidates are required to have pj > 25 GeV and either I 77 I < 1.37 or 1.52 < \r]\ < lAl to 
exclude the transition between the barrel and endcap calorimeters. Muon candidates are required to have 
Pt > 25 GeV and It/I < 2.5. Nonprompt leptons and nonleptonic particles may be reconstructed as leptons 
and satisfy the selection criteria, giving rise to nonprompt and fake lepton backgrounds. In the case of 
electrons, these include contributions from semileptonic decays of b- and c-quarks, photon conversions, 
and jets with a large electromagnetic energy fraction. Nonprompt or fake muons can originate from 
semileptonic decays of b- and c-quarks, from charged-hadron decays in the tracking volume or in had¬ 
ronic showers, or from punchthrough particles emerging from high-energy hadronic showers [54]. The 
nonprompt and fake lepton backgrounds are reduced by requiring the lepton candidates to be isolated 
from other energy deposits or high-px tracks. The tracks used in the isolation calculation are required to 
originate from the primary interaction vertex and have px > 1 GeV. For electrons, an //-dependent limit 
is placed on the amount of energy measured in the calorimeter within a A/? = 0.2 cone around the can¬ 
didate which is neither from the electron candidate itself nor from additional pp interactions. A similar 
requirement is placed on the scalar sum of the px of tracks within a A/? = 0.3 cone around the track of 
the electron candidate. Each requirement has an average efficiency of 90% for electrons from Z ^ ee. 
Muons are required to have a px-dependent track isolation [62], requiring the scalar sum of the px from 
tracks with A/? < 10 GeV//r^ to be less than 0.05 • p^, where //^ is the px of the candidate muon track. 
This requirement has an efficiency of approximafely 95% for muons from W pv. 

In this analysis, r leptons are not explicitly reconstructed. Because of the high px threshold, only a small 
fraction of r leptons decaying leptonically are reconstructed as electrons or muons, while the majority of 
T leptons decaying hadronically are reconstructed as jets. 


Jet requirements 

Events must contain at least four jets with px > 25 GeV and |?/| < 2.5 reconstructed using the anti-kf 
algorithm [63] with a radius parameter R = 0.4. The jets are constructed from calibrated topological 
clusters built from energy deposits in the calorimeters, and they are calibrated to the hadronic scale [64]. 
Prior to jet finding, a local cluster calibration scheme [65] is applied to correct the topological cluster 
energies for the effects of calorimeter noncompensation, dead material and out-of-cluster leakage. The 
jet energy scale was determined using information from test-beam data, EHC collision data, and simula¬ 
tion [64, 66 ]. All jets are required to have at least two tracks that originate from the primary interaction 
vertex. In order to suppress jets that do not originate from the primary vertex, jets with px < 50 GeV and 
\r]\ < 2.4 are required to have a jet vertex fraction (JVF) above 0.5, where JVF is the ratio of the summed 
scalar px of tracks originating from the primary vertex to that of all tracks associated with the jet. An 
overlap removal procedure [61] is applied to remove jets that were already identified as electrons. 

To identify jets originating from the hadronization of a fr-quark (‘fr-tagging’), a continuous discriminant 
is produced by an algorithm using multivariate techniques [67] to combine information from the impact 
parameter of displaced tracks and topological properties of secondary and tertiary vertices reconstructed 
within the jet. The efficiency for a jef confaining a fr-hadron fo be fr-fagged is 70%, while fhe lighf-jef (c- 
jef) efficiency is less than 1 % ( 20 %) as determined in simulated tt events, where fight jets are jets initiated 
by a U-, d-, 5 -quark, or gluon. If any jets are fr-tagged, the event is rejected. 
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Missing transverse momentum requirements 


The is constructed from the vector sum of calibrated energy deposits in the calorimeter and re¬ 
constructed muons [68]. Events that do not contain a leptonically decaying W boson are suppressed by 
requiring > 30 GeV and (E™®® -i- ) > 60 GeV, where is the transverse mass of the W boson 

defined as trij = - cos A(p), where A0 is the azimuthal angle between the charged-lepton 

transverse momentum vector and 

Table 1: Observed and expected event yields after preselection. The quoted uncertainties include both the statistical 
and systematic contributions, with the latter being dominant. The sources of systematic uncertainty are discussed 
in Sec. 6. 



Electron 

Muon 

VT-i-jets 

110000 So 

145000 +20000 

iH-juuu _28000 

Z-i-jets 

28000 

15200 "^2200 

1 jzuu 

tt 

19200 

ivz,uu _2900 

23700 

Diboson 

3900 ±1900 

4400 ±2200 

Single top 

3750 

-560 

4500 !®o 

Multijet 

22000 ±8800 

13300 ±5300 

Total background 

183000!™ 

206000 +2™oo 

zuuuuu 34000 

Signal {niQ - 700 GeV) 

79 ±10 

79 ±10 

Data 

182075 

208641 


4.2 Final selection 

With the final-state objects identified, additional kinematic requirements are applied to exploit the distinct 
features of the signal, assumed to be heavier than the previously excluded mass of 350 GeV. The W boson 
and light quark originating from Q Wq would be very energetic and have a large angular separation 
due to the large mass of the Q quark. On the other hand, the decay products of the W boson would tend to 
have a small separation due to the W boson’s boost. By selecting only events that are consistent with these 
properties, the VT-i-jets background yield is reduced by orders of magnitude. To facilitate the discussion 
of the kinematic selection, the following objects are defined: 

• Bhad is the candidate for a W boson in the decay Q Wq qqq', 

• Vkiep is the candidate for a W boson in the decay Q —> Wq —> ivq\ 

• ( 7 -jet is a candidate for the jet originating from the qinQ ^ Wq (i.e., from the decay of the heavy 
quark, not the hadronically decaying W boson). There are two ^-jets per event, so qi (^ 2 ) is used to 
denote the one with higher (lower) pT- 
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The Whad candidate is defined as a dijet system with pj > 200 GeV, angular separation AR < 1.0, and 
an invariant mass in the range 65-100 GeV. All possible jet combinations are considered and, if multiple 
pairs satisfy the above requirements, the pair with the mass closest to that of the W boson [69] is chosen. 
If no VThad candidate is found, the event is removed. This requirement results in 94% background rejection 
while keeping 53% of the signal if mq = 700 GeV. The mass distribution of the VThad candidates prior to 
the mass requirement is shown in Fig. 2. 




mjj [GeV] mjj [GeV] 


Figure 2: Invariant mass distributions for dijet systems with px > 200 GeV, angular separation AR < 1.0, and an 
invariant mass closest to the W boson mass, for the (left) electron and (right) muon channels after the preselection 
requirements. Vkhad candidates correspond to the mass range 65-100 GeV. The highest bin includes all events with 
ifijj > 240 GeV. The shaded band shows the total uncertainty on the background prediction. The lower panel shows 
the significance of the difference between data and expectation in units of normal standard deviations as explained 
in Ref. [70]. 


The VFiep candidate is reconstructed using the lepton and which is taken to be the neutrino pj. 

The longitudinal momentum of the neutrino is determined up to a two-fold ambiguity by requiring the 
invariant mass of the lepton-neutrino system to equal the mass of the W boson. When no real solution 
exists, the neutrino pseudorapidity is set equal to that of the lepton because the decay products of the 
W boson tend to be nearly collinear for the kinematic regime of interest. In simulated samples, the rate 
of events with no real solution is approximately 35%. Signal events are expected to have energetic W 
bosons, so the VTiep candidate is required to have pr > 125 GeV. 

The W candidates (Whad and Wiep) are then each paired with a different one of the remaining jets to 
create the two heavy quark candidates, Q1 and Q2. This step involves testing all possible pairings of 
^-jet candidates with the VFhad and VTiep candidates. In addition, the Wiep candidate may have two real 
solutions for the longitudinal momentum of the neutrino. Among the possible combinations of neutrino 
momentum solutions and Wq pairings, the one yielding the smallest absolute difference between the two 
reconstructed heavy quark masses is chosen. In simulated samples, the rate of correct Wq pairing is 
approximately 40% (48%) for a signal of mass 400 GeV (800 GeV). Once the heavy quark candidates are 
determined, the ^-jets are required to have pxCvi) >160 GeV and pt( 92 ) > 120 GeV and the difference 
between the reconstructed heavy quark masses must be less than 120 GeV. 






















With candidate objects identified for two heavy quarks and their decay products, the following additional 
kinematic criteria are applied. Each event must have //j > 1100 GeV, where //t is the scalar sum of the 
transverse momenta of the lepton, neutrino, Whad> and the two i^-Jets. The angular separation between the 
lepton and neutrino candidates must satisfy A7?(^, v) < 1.4. The heavy quarks would tend to be central 
and back-to-back, so the angular separation between the two reconstructed heavy quarks is required to 
satisfy 2.0 < AR{Ql, Q2) < 4.2. 

To further capitalize on the presence of a hadronically decaying W boson with high pj in the signal, the 
analysis makes use of a splitting variable [30] defined as: 

mm{pTji,pjj2f X AR(juj2f 

yn = - 5 -; 

m , 

where j\ and 72 are fhe fwo jefs from the Whad candidate and mjij 2 is the mass of the Whad candidate. 
The two jets from a hadronically decaying W boson tend to have roughly equal energies, while dijets 
from QCD processes are likely be asymmetric in energy. Furthermore, jets from W ^ qq tend to have a 
larger opening angle due to the mass of the W boson. Dividing by the dijet mass provides discrimination 
for cases in which a QCD dijet system happens to have large values for min(pT/i> and AR{j\, 72 )^, 
as such systems are likely to have a very large mass. The background is reduced by a factor of 3.2 by 
requiring yu > 0.25. This requirement has a signal efficiency of approximately 50%, although the precise 
value depends on the mass of the heavy quark. The modeling of the splitting variable with a large number 
of events can be seen in the top histograms of Fig. 3, which depict the yi 2 distribution after preselection 
and requiring a IThad candidate. The bottom of Fig. 3 contains the distributions of yi 2 immediately before 
the yi 2 > 0.25 requirement is applied. 

The final selection criteria are motivated by the fact that the decay products from the Q quark are well 
separated. The final requirements are A7?(lThad, <?i) > 1-0, A7?(lThad, <? 2 ) > 1-0, A/?(lTiep,^i) > 1.0, and 
A/?(VTiep, < 72 ) > TO. Table 2 presents a summary of the expected and observed numbers of events after the 
final selection, for which the signal (background) efhciency compared to preselection is approximately 
8 % (0.004%). The small contributions from tt, Z-i-jets, dibosons, single-top, and multijet events are 
combined into a single background source referred to as ‘non-lT-i-jets’. Uncertainties on the yields include 
the uncertainty due to the size of the signal and background samples and the cumulative effect of the 
systematic uncertainty described in Sec. 6 . 

The final discriminant variable used in this search is nij-eco, the reconstructed heavy quark mass built from 
the VThad candidate and the paired ^-jet candidate. 


5 Background modeling 

5.1 Correction to W-i-jets modeling 

It is observed after applying the preselection criteria that the simulated VT-i-jets sample does not accurately 
model the px spectrum of the leptonically decaying W boson candidate. This mismodeling leads to an 
overestimation of the IT-i-jets yields in the high-momentum tails of the F™®®, lepton px. jet Pt, and 
Hj distributions. The dominant background for this analysis is W-i-jets events in which the transverse 
momentum of the W boson, px(lT), is high. Therefore, it is important to have accurate predictions for 
both the overall normalization and the px(lT) distribution. This section describes the procedure used to 
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Figure 3: The splitting variable yn between the decay products of the hadronic W boson (top) after preselection 
and requiring a hadronic W boson candidate and (bottom) immediately before applying the requirement yi 2 > 0.25, 
for the (left) electron and (right) muon channels. The shaded band shows the total uncertainty on the background 
prediction. The lower panel shows the significance of the difference between data and expectation in units of normal 
standard deviations. 
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Table 2: Expected yields for the backgrounds and the VLQ signal with mg — 700 GeV, along with the observed 
number of data events, after applying all selection criteria. The uncertainties on the predicted yields correspond to 
the statistical uncertainty due to finite sample size and the systematic uncertainty, respectively. 



Elecfron 

Muon 

VF-i-jels 

Non-VF-i-jefs 

5.6±1.51};5 
1.2 ±0.5 

6.0 ±1.0 +2-2 

1.2 ±0.4 

Tofal background 

6 . 8 ± 1 . 6 : 2.4 

7.2±1.1!2J 

Signal {niQ = 700 GeV) 

7.0 ±0.6 +}] 

6.9 ±0.6 +} ° 

Dafa 

9 

11 


derive the VT+jets and multijet normalizations and a reweighting function to correct the vector-boson pj 
in the VT-i-jets and Z-i-jets^ simulated samples. All steps in the procedure rely on fits to the data using the 
Pt distribution of the VTiep candidate, after applying the preselection requirements. 

First, the normalizations for the VF-i-jets templates and the multijet templates are fit to the data, with 
all other background processes fixed af fheir expecfafions. The difference befween fhe observed and 
predicfed number of IF-i-jefs evenfs is assumed fo be due fo fhe cross secfion differing from ifs predicfed 
value, so fhe elecfron and muon channels are fil simulfaneously fo defermine a single fF-i-jefs scale faclor 
of 0.82. A correction for is fhen derived fhaf minimizes fhe befween dafa and simulation 

for fhe p!p®‘^°(VFiep) disfribufion, where is fhe pj of fhe generated vecfor boson in fhe VF-i-jels or 

Z-i-jefs sample jusf before ifs lepfonic decay. The reweighfing is approximafely unify for low 
decreasing fo 0.86 for = 200 GeV and 0.65 for = 500 GeV. Finally, fhe normalizations 

for fhe mulfijef samples and fhe corrected fF-i-jefs samples are til fo fhe dafa, wifh all ofher background 
processes fixed. The fif is done simulfaneously in fhe elecfron and muon channels. Figure 4 shows fhe 
^reco(jyiep) distribufion in fhe elecfron and muon channels afler applying fhe p!j7’*(V) correcfion and scale 
facfors. The corresponding disfribufions for //t and are shown in Figs. 5 and 6. The uncerfainfies 
on fhe normalizafions and reweighfing are described in Sec. 6. 

5.2 Validation regions 

The following validafion regions are used fo verify fhe modeling of fhe background processes: 

• VRl: preselection, plus one fFhad and ntreco < 350 GeV; 

• VR2: preselecfion, plus one W^ad and Hj < 800 GeV ; 

• VR3: final selection, buf wifh fhe requiremenfs on priqi), PTiqi), and Hj changed fo pjiqi) < 
160 GeV, pt(^ 2 ) < 80 GeV, and Hj < 800 GeV and wifh no requiremenfs on A/?(VFhad>^i)> 
AR(iVhad,^ 2 ), AR(iViep,^i), and AR(iViep, ^ 2 )- 

The validafion regions are orthogonal fo fhe signal region, buf are nevertheless useful for checking wifh a 
larger number of evenfs fhaf fhe background normalization and kinematics are well-modeled. The expec- 
fed signal confribufion in each validafion region is approximafely fhe size of fhe background uncerfainfy 


^ The correction is primarily motivated by the dominant Vf+jets background, but is also applied to the Z+jets background. 
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Figure 4; The transverse momentum distributions for the leptonic W boson candidate after the correction 

and hnal normalization for the (left) electron and (right) muon channels. All overflows are included in the right¬ 
most bin. The shaded band shows the total uncertainty on the background prediction. The lower panel shows the 
signihcance of the difference between data and expectation in units of normal standard deviations. 
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Figure 5: Distributions of the scalar sum of the transverse momenta of the lepton, neutrino, Whad, and the two 
q'-jets (Hj) after the pl^“*(y) correction and hnal normalization for the (left) electron and (right) muon channels. 
All overhows are included in the rightmost bin. The shaded band shows the total uncertainty on the background 
prediction. The lower panel shows the signihcance of the difference between data and expectation in units of normal 
standard deviations. 


for signal masses around the previous limit (350 GeV) and decreases very rapidly for higher masses. VR2 
and VR3 are used to validate the modeling of the final discriminant, nij-^co, as shown in Figs. 7 and 8. 
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Figure 6; Distributions of the missing transverse momentum after the correction and hnal normalization 

for the (left) electron and (right) muon channels. All overflows are included in the rightmost bin. The shaded band 
shows the total uncertainty on the background prediction. The lower panel shows the signihcance of the difference 
between data and expectation in units of normal standard deviations. 


VRl is used to verify the modeling of variables other than nireco, such as the distribution and the pj 
spectra for the individual objects. 




mreco [GeV] [GeV] 

Figure 7: Comparison between data and simulation for the distribution of nireco, the reconstructed heavy quark 
mass built from the Whad candidate and the paired q'-jet candidate, in validation region VR2 for the (left) electron 
and (right) muon channels. All overflows are included in the rightmost bin. The shaded band shows the total 
uncertainty on the background prediction. The lower panel shows the signihcance of the difference between data 
and expectation in units of normal standard deviations. 
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Figure 8; Comparison between data and simulation for the mreco distribution in validation region VR3 for the (left) 
electron and (right) muon channels. The shaded band shows the total uncertainty on the background prediction. 
The lower panel shows the significance of the difference between data and expectation in units of normal standard 
deviations. 


6 Systematic uncertainties 


The uncertainties considered in this analysis can affect the normalization of signal and background and the 
shape of the final discriminant, nxreco- Each source of systematic uncertainty is assumed to be 100% cor¬ 
related across all samples, but the different sources are treated as uncorrelated with one another. Table 3 
shows the impact of the dominant uncertainties on the normalization of the background processes and a 
signal sample with a mass of 700 GeV. 

Table 3; Table summarizing the overall normalization changes (expressed in %) in signal and background yields 
for the dominant systematic uncertainties considered. The selection presented here is the combination of the e-i-jets 
and yu-Hjets channels after the final selection. 


Signal {niQ - 700 GeV) Non-VF-i-jets IV-i-jets 


Luminosity 

-t2.8/-2.8 

-t2.8/-2.8 

±2.8/-2.8 

Normalization 


±15 

±2.7/-4.4 

Lepton identification 

+1.6/-1.6 

±1.5/-1.5 

±1.4/-1.4 

Jet energy resolution 

+0.6/-0.6 

±12/-12 

±8.7/-8.7 

Jet energy scale 

+6.1/-4.3 

±33/-34 

±14/-18 

f7-tagging 

+0.11-0.1 

±5.1/-5.3 

±0.3/-0.3 

c-tagging 

+1.5/-1.5 

±1.5/-1.5 

±1.2/-1.2 

Light-jet tagging 

+1.0/-1.0 

+0.91-0.9 

±1.0/-1.0 

^ituth(y) re-weighting 



+5.1I-A.1 
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6.1 Normalization uncertainties 


Uncertainties affecting only the normalization include those on the integrated luminosity (±2.8%) and the 
cross sections for various background processes. The uncertainty on the integrated luminosity is derived 
following the same methodology as that detailed in Ref. [71]. This uncertainty is applied to all simulated 
signal and background processes. 

After the final selection, the non-lT±jets background has a total normalization uncertainty of 15%. The 
predicted contribution from the multijet background is negligible compared to the uncertainty on the 
non-Vk±jets background, so it is neglected. 

As described in Sec. 5.1, the normalization of the lT±jets background is determined from a fit to data 
using both lepton channels. The uncertainty on the lT±jets normalization is determined by comparing 
that result to the normalization one would obtain if only the electron or only the muon channel were used. 
This is motivated by the fact that it is not known whether the normalization from the electron channel 
or the muon channel (or something in between) is correct and this procedure leads to an uncertainty of 
±2.7/-4.4%. The statistical uncertainty from the fit is negligible (0.03%) in comparison. 

The rest of the systematic uncertainties can modify both the normalization and shape of the mreco distri¬ 
bution. 


6.2 Shape uncertainties 

Uncertainties on the trigger, reconstruction, and isolation efficiencies for the selected lepton are estim¬ 
ated using Z ^ ee and Z ^ fi/u events [72, 73]. In addition, high jet-multiplicity Z ^ events are 
studied, from which extra uncertainties on the isolation efficiency are assigned to account for the differ¬ 
ence between Z boson and tt events. Uncertainties on the reconstruction and the energy scale and 
resolution of the leptons were also considered; however, these have a very small impact on the results. 

The jet energy resolution is measured by studying dijet events in data and simulation. The simulation is 
found to agree with data to better than 10% [74]. The differences in resolutions between data and simu¬ 
lations are used to determine the relative systematic uncertainty. The uncertainty on the jet energy scale 
is evaluated by repeating the analysis with the jet energy scale shifted by ±10" [64, 66]. The jet recon¬ 
struction efficiency is estimafed using track-based jets and is well described in simulation. To account 
for differences in the efficiency for reconstructing jets in simulated events compared to collider data, the 
efficiencies are measured in both sampies and the difference is taken as the uncertainty. The uncertainty 
due to the JVF requirement is evaiuated by comparing the signai and background distributions with the 
JVF cut shifted up and down by i0%, a variation that spans the difference observed befween data and 
simuiation in this quantity. The (^-tagging efficiency for (j-jefs, as well as c-jefs and light jets, is derived 
in data and a simulated tt sample, parameterized as a function of pj and rj [67, 75]. The corresponding 
(in)efiiciencies in the simulated samples are corrected to match those in data and the uncertainties from 
the calibration are propagated through the analysis. 

The VT±jets sample is assigned a pI^'“*(U)-dependent shape uncertainty due to the correction described in 
Sec. 5.1. Four sources of uncertainty on the p™^(U) correction are considered: (1) the statistical uncer¬ 
tainty on the parameters of the reweighting function; (2) the difference between the nominaf reweighting 
function and alternative corrections obtained by only considering the electron or muon channel; (3) the 
dependence of the fit on the choice of bin width when deriving the reweighting function; (4) the difference 


i5 


between alternative parameterizations of the reweighting function. A closure test is performed to verify 
that any residual differences between data and prediction are well within the assigned uncertainty. 


7 Results 


The final mreco disfribution for fhe combined elecfron and muon channels is shown in Fig. 9 for fhree 
signal scenarios: tuq = 600 GeV, iuq - 700 GeV, and tuq - 800 GeV. The observed disfribution shows a 
slight excess over the SM expectation, but the excess is broader than expected for signal and is consistent 
with the background-only prediction at the level of two standard deviations. Therefore, the analysis 
proceeds to setting limits on the signal hypothesis. 

The nireco distribution for the combined electron and muon channels after the final selecfion (Fig. 9) 
is used fo derive 95% conhdence level (CL) limifs on fhe QQ producfion cross secfion using fhe CLs 
mefhod [76, 77]. 
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Figure 9: Distribution of the final discriminant, the mass of the hadronically decaying heavy quark candidate, with 
the bin widths chosen as for the statistical analysis. The expected contribution from signal with BR(2 ^ Wq)-l is 
shown stacked on top of the total background prediction for three mass scenarios, otq = 600,700, and 800 GeV. The 
backgrounds are stacked with the largest on top. The shaded band shows the total uncertainty on the background 
prediction. 


Limits on the pair production of new chiral quarks are evaluated by setting BR(2 ^ Wq) = 1. Fig¬ 
ure 10 shows the observed and expected limits on a heavy chiral quark as a function of niQ, compared 
to the theoretical prediction [33-38]. The total uncertainty on the theoretical cross section includes the 
contributions from the scale variations and PDF uncertainties. Using the central value of the theoretical 
cross section, the observed (expected) 95% CL limit on the mass of a new chiral quark is mq > 690 GeV 
(780 GeV). This represents the most stringent limit to date on the mass of a new quark decaying exclus¬ 
ively into a W boson and a light quark (u, d, s). This limit is also applicable to the production of pairs of 
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Figure 10; The expected (black dashed line) and observed (black solid line) 95% CL upper limits on the cross 
section as a function of itiq when setting BR((2 —» Wq) - 1, which would be the case for a new chiral quark. The 
green and yellow shaded bands indicate the +lcr and +2cr intervals on the distribution of expected results for the 
chiral model if no signal exists. The blue line is the theoretical prediction for the total cross section [i.e., assuming 
BR((2 —> Wq) = 1] with its uncertainties. 


down-type vectorlike quarks with electric charge of -4/3 which each decay into a W boson and a light 
quark (d, s). 

Next, the VLQ signal samples are used to set limits on the mass of a heavy quark that decays to a light 
quark (u, d, s) and either a W, Z, or H boson. The results are given as a function of the branching 
ratios BR(2 ^ Wq) versus BR(2 ^ Hq), with the branching ratio to Zq fixed by the requirement 
BR(2 —> Zq) = 1 - BR(2 ^ Wq) - BR(2 ^ Hq). The analysis loses sensitivity at low masses due to 
the tight selection requirements optimized for the decay Q —> Wq, so the results are presented as the 
upper and lower bounds on the mass range that is excluded at 95% CL. The expected limits on mg as a 
function of the branching ratios are shown in Fig. 11 and the observed limits are shown in Fig. 12. For 
example, for the branching ratios BR(2 ^ Wq) - 0.6, BR(2 ^ Zq) - 0.4, and BR(2 ^ Hq) - 0, 
the results exclude VLQs with a mass from 370 GeV to 610 GeV. If no signal is present, the median 
expected exclusion range for that set of branching ratios is 340 GeV to 690 GeV. For BR(2—> IFi^) = 0.6, 
BR(2 ^ Zq) = 0.3, and BR(2 ^ Hq) = 0.1, the results exclude VLQs with a mass from 370 GeV to 
470 GeV, while the median expected exclusion range is 340 GeV to 680 GeV. For BR(2 ^ Wq) = 0.5, 
BR(2 —> Zq) = 0.25, and BR(2 ^ Hq) - 0.25, the results exclude VLQs with a mass from 390 GeV to 
410 GeV, while the median expected exclusion range is 360 GeV to 650 GeV. 
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Figure 11: The (left) upper and (right) lower bounds on the range of heavy quark masses expected to be exclude 
at 95% CL, as a function of the branching ratio of the heavy quark to Wq versus Hq, with the branching ratio t 
Zq fixed by the requirement BR(2 —> Zq) = 1 - BR(2 —> Wq) - BR(2 —> Hq). The region above the diagonal i 
forbidden by unitarity. 
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Figure 12: The (left) upper and (right) lower bounds on the range of heavy quark masses observed to be excluded 
at 95% CL, as a function of the branching ratio of the heavy quark to Wq versus Hq, with the branching ratio to 
Zq fixed by the requirement BR(2 Zq) = 1 - BR(Q Wq) - BR(2 ^ Hq). The region above the diagonal is 
forbidden by unitarity. 


8 Conclusion 

A search for new heavy quarks that decay to a VT boson and a light quark (u, d, s) is performed with a 
dataset corresponding to 20.3 fb“^ that was collected by the ATLAS detector at the LHC in pp collisions 
at = 8 TeV. No significant evidence of a heavy quark signal is observed when selecting events with 
one charged lepton, large and four non-fj-tagged jets. The results exclude new heavy chiral quarks 
with masses up to 690 GeV at 95% CL. This search is also interpreted in the context of vectorlike quark 
models, for which the new heavy quark can decay to a light quark and either a W, Z, or H boson. New 
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exclusions limits are presented in the two-dimensional plane of BR(2 ^ Wq) versus BR(2 ^ Hq). 
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